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Abstract 

We investigate on the lattice the metastable confined phase above Tc in SU (N) gauge theories, 
for = 8, 10, and 12. In particular we focus on the decrease with the temperature of the mass of 
the lightest state that couples to Polyakov loops. We find that at T = Tc the corresponding effective 
string tension o"efr(T) is approximately half its value at T = 0, and that as we increase T beyond 
Tc, while remaining in the confined phase, cJefr(T) continues to decrease. We extrapolate o"efr(T) 
to even higher temperatures, and interpret the temperature where it vanishes as the Hagedorn 
temperature Th- For SU{12) we find that Th/Tc = 1.116(9), when we use the exponent of the 
three-dimensional XY model for the extrapolation, which seems to be slightly preferred over a 
mean-field exponent by our data. 

PACS numbers: 11.15.Ha,11.15.Pg,11.25.Tq 
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I. INTRODUCTION 



In the large-iV limit the confined phase of the SU{N) gauge theory becomes weakly 
interacting and relatively simpler [l|, |^. Moreover, some of its features can be described 
by a low energy effective string theory (for example see j^). The link between gauge theo- 
ries and string theories was strengthened with the conjectured AdS/CFT dualities between 
supersymmetric Yang-Mills theories in the large-A^ limit and gravity models (for a review 
see 0). An interesting element of systems with stringy properties is the occurrence of an 
ultimate temperature, above which the free string description is unsuitable. In string models 
for QCD (or pure gauge theory), this temperature is naturally identified with a 'Hagedorn' 
deconfinement transition. 

The earliest evidence for the existence of a finite temperature phase transition in hadronic 
physics was obtained by Hagedorn in There, he assumed that hadronic states with mass 
m —>■ oo, referred to as "fireballs", are compound systems of other fireballs. This can be 
stated mathematically by a self-consistent relation that the density of states p{E) must obey, 
and that results in an exponential growth of p with the energy E. An immediate consequence 
is that there exists a temperature T^, above which the partition function diverges: as 
T — > T^, the Boltzmann suppression of states oc exp{—E/T) is overwhelmed by the density 
of states, p{E) oc exp(+ci?). At this T an increase of the system's energy is met with an 
increase in the number of particles. In the bootstrap dynamical framework, which predated 
the discovery of QCD and the idea that hadrons are confined bound states of quarks and 
gluons, Th represents an ultimate temperature beyond which matter cannot be heated. In 
the more appropriate language of QCD, once we are close enough to Th for these fireballs 
to be densely packed, the underlying quark and gluon degrees become liberated from their 
'hadronic bags' and we expect to see a deconfinement transition 

A simple, intuitive, and general analogue of the above argument 0, 0] context of 

any linearly confining theory is as follows. In such a theory the energy of a string of length 
/ between two distant static sources a distance r < / apart, obeys 

E(/;r)~a/, (1.1) 

where a is the confining string tension and we neglect subleading terms. For a string (or for 
a long enough fiux tube) it is easy to see that once / ^ r the number of different states of 
the string grows exponentially with /, n{l) oc exp(c/), up to power factors, with c determined 
by the dynamics and dimensionality. On the other hand the probability of such a string 
with I ^ r is suppressed by a Boltzmann factor, oc exp{—E{l)/T). The total probability of 
such a string is given by the free energy that combines the two factors to define an effective 
string tension (Teff(T): 

oc exp [(c - a/T) ■ I] = exp {-a,s{T)l/T) . (1.2) 

It is thus clear that as T Th = cr/c, arbitrarily long loops will be thermally excited and 

the effective string tension vanishes, aes — " 0. Because of this vanishing energy, and other 
features of the free energy , one would naturally expect this deconfining phase transition 
to be second order. 

We can give a parallel argument for the partition function itself. Consider an SU{N) 
gauge theory in which the confining energy eigenstates are composed of glueballs. A model 
for glueballs is to construct them from closed loops of (fundamental) fiux. For the lightest 
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glueballs this loop will be small, and given that the width of the flux tube is 0(1/a/(t), it 
does not make much sense to think of a distinct closed loop of flux. However in the sector of 
highly excited glueballs, the loop will be very long and the presence of such states becomes 
compelling. For a state composed of a flux loop of length I, the energy is ~ al. What is 
the number density, ndE), of such states? Such highly excited states have large quantum 
numbers and in this limit a classical counting of states is justified. Thus the number of these 
states equals the number of closed loops of length I, i.e. nc^E) oc exp(c/) up to subleading 
factors, with the scale c identical to the one in Eq. ()1.2j) . Thus the partition function will 
have a Hagedorn divergence as T — > Th where Th is identical to the string condensation 
temperature of the previous paragraph. 

Since we expect c ~ ^Ja, we expect the Hagedorn transition to occur at = a /c ^ ^Ja. 
On the other hand the lightest glueball (a scalar) satisfies mc — 3.5^^(7 while the next lightest 
glueball (a tensor) satisfies mc — 4.8^0". Thus one has the somewhat counterintuitive 
picture that as T — Th the lightest glueballs are not thermally excited. Rather it is the 
highly excited glueballs, whose density of states grows exponentially with their mass, that 
drive the transition. 

In the above scenarios, as T — the vacuum becomes increasingly densely packed with 
the thermally excited states. These will at some point start to interact and the idealised 
arguments we use necessarily break down as T approaches Tfj. We note that as — oo for 
SU{N) gauge theories (or QCD), interactions between colour singlet states vanish. Thus it 
is in this limit that the argument for a Hagedorn transition becomes most compelling. 

The vanishing of aes, and the divergence of the associated correlation length, suggests 
that this Hagedorn transition is second order. To what universality class might it belong? 
The high temperature phase has a nonzero vacuum expectation value for the complex valued 
Polyakov lines that wind around the Euclidean temporal torus. This spontaneously breaks 
the global symmetry of the theory. Using universality arguments one can then predict 
the critical exponents of the transition In particular, in three spatial dimensions it 
belongs to the universality classes of the three-dimensional Ising and XY models for N = 2 
and > 4 (for = 3, there is no known universality class) [lO|. Finally for N = oo there 
are studies that predict mean- field behaviour for the correlation length This can be 

understood as a suppression of the critical region by powers of (see for example the 
discussion in [l^], and references therein). The study in 113] also gives mean- field scaling, 
however only because infrared divergences are not seen at the small volume discussed there 

a. 

The above discussion has so far ignored the contribution to the partition function that 
comes from nonconfined energy eigenstates containing a finite density of gluons. While such 
states will be irrelevant at low T, the fact that their entropy grows as A^^ while the entropy 
of confined states is at most weakly dependent on A^, means that for large enough A^ there 
must be some T where their free energy will decrease below that of the confined sector of 
states. At this point there will be a phase transition to a (perhaps strongly interacting) 
'gluon plasma', which one would naturally expect to be first order. Indeed it turns out to be 
the case that in four dimensions SU{N) gauge theories go through a first order deconfining 
transition for A^ > 3 (See also the latent heat calculation at large- A^ on a symmetric 
lattice For SU(2) the transition is second order, but it is not clear if it is a Hagedorn 

transition. The fact that it is at small rather than large A^ makes the case weaker. As 
does the fact that the SU{2) value of T^j \/o lies on the curve that interpolates through the 
A^ > 3 values On the other hand the value of Td^Jo does coincide with the string 
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condensation temperature of the simple Nambu-Goto string theory (see below). 

The fact that for > 3 the first order deconfining transition occurs for < Th, would 
appear to render the Hagedorn transition inaccessible. However the deconfining transition 
is strongly first order at larger A^, and so one can try to use its metastabihty to carry out 
calculations in the confining phase for T > Tc. If Th is not far away, one can then hope 
to calculate aes{T) over a range of T where it decreases sufficiently that an extrapolation 
to creff(Tff) = can be attempted. In the range of N accessible to us, the interface tension 
between confined and deconfined phases increases with faster than the latent heat, and 
this makes the metastabihty region larger [17] as N grows. Thus such a strategy has some 
chance of success, and it is what we shall attempt in this paper. 

Our strategy is therefore to begin deep enough in the confined phase and then to increase 
the temperature to temperatures T > Tc, calculating the decrease in aes(T), and extrapo- 
lating to (Teff = 0. We interpret the result of the extrapolation as the Hagedorn temperature, 
Th- Nevertheless, since we work with finite values of N and volume V, tunneling probably 
occurs somewhere below Th- These tunneling effects and the fact that as ciefr decreases finite 
volume effects become important, can make an apriori fit for the functional form of aes{T) 
unreliable. As a result we first perform fits where we fix the functional behaviour to be 

m = oTeff/T --J« {Th/T, - T/T,Y . (1.3) 

where u = 0.6715(3) corresponding to 3D XY or u = 0.5 corresponding to mean field. 
The reason for these two choices is motivated by two conceivable ways in which the low 
energy effective loop potential can behave (see below). In addition we also perform fits 
where the exponent z/ is a free parameter, constrained to the range [0, 1], and find it to be 
especially useful for SU{12)- The coefficient A is fitted as well. An additional outcome of 
this work is to confirm that at T = Tc the mass of the timelike fiux loop that couples to 
Polyakov loops is far from zero at large-A^, which confirms that the transition is strongly 
first-order. To make this point clear we will present figures of the effective string tension 
in units of the zero temperature string tension ex. As a function of T, it should have the 
following behavior. 

a^s/a T^J« A ■ T/T, ■ {Th/T, - T/T^ , (1.4) 

If we imagine the effective potential for an order parameter such as the Polyakov loop, 
Vee{lp), then aX T = T^ this will possess degenerate minima corresponding to the confined 
and deconfined phases. These will be separated by a barrier whose height is expected to be 
0(A^2) at large A^. As we increase T the confined minimum rises relative to the deconfined 
one(s). As T — > Th we expect the second derivative at the confining minimum to go to zero, 
corresponding to (Jefr 0. The simplest possibility is that at this T the confining minimum 
completely disappears, i.e. that this corresponds to a spinodal point of the potential. In a 
string model of glueballs where the glueball is composed of two 'constituent' gluons joined by 
an adjoint fiux tube, string condensation would correspond to the explicit release of a gluon 
plasma simultaneously throughout space and the identification of the Hagedorn transition 
with a spinodal point would be compelling. This is less clear in the closed flux loop model. 
(Whether the increasing stability of the adjoint string as A^ — > cxo allows us to use the adjoint 
string model for the highly excited states relevant near Th also needs consideration.) If this 
Hagedorn transition-spinodal point identification is indeed correct, then in the vicinity of 
Ip ~ 0, the effective loop potential looks like that of a Gaussian model, and v = 0.5 is 
expected. A model for this behaviour is in where because of the infinitesimal volume. 
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any infra-red divergences are excluded, and one has mean-field scaling close to Th-, which 
also implies v = 0.5. 

In principle it seems quite possible that does not coincide with the spinodal transition 
temperature, Tg. The temperature Th is a natural concept if one has a good description of 
the confined phase as an effective string theory. The latter will have 0{1/N) interactions, 
and thus to leading order, lead to a Hagedorn behaviour at a certain temperature Th- On 
the contrary, Tg most probably encodes information on the gluonic deconfined phase, which 
might not be contained in the string theory. Without any other information in the spirit of 
the calculations 13], that identifies Th = Tg, these two temperatures may be different. 

If Th < Tg then our method for identifying Th remains valid. In that case one may write 
down a Landau- Ginzburg theory that has a second order Hagedorn transition at Th, from the 
confined vacuum, C, to a deconfined one, Di. This happens at Th > Tc, where both C, Di are 
metastable, and another deconfined vacuum D2 is stable. At large- iV, the metastability may 
be strong enough such that this embedded second order phase transition happens without 
being sensitive to tunnelings into the real vacuum. In this case, the fixed point that controls 
the critical behavior of the transition is that of the 3D XY model. Nonetheless, to see this 
nontrivial critical behavior one must be very close to Th- This happens because interactions 
between the Polyakov loops are 0(l/iV^) suppressed, and taking the N = 00 limit before 
T Th results in a Gaussian model and to u = 0.5. To see the scaling behavior of the 
3DXY model, one needs to take T — * Th together with — >■ 00. In other words, the critical 
region of this second order transition has a width that shrinks with increasing N jl^ . To see 
this we can examine the renormalization of the A|/p|^ interaction in the Landau-Ginzburg- 
Wilson effective action of the transition. As usual, A gets renormalized by infra-red (IR) 
modes, and in three spatial dimensions it gets a contribution of 0(A^/m). Since A ~ 1/A^^, 
this becomes significant only if m ~ l/N"^. Using ~ {Th — T), this means that the IR 
modes drive the system to the 3DXY universality class only if {Th — T)< Outside of 

this regime the correlation length has the MF scaling u = 0.5. 

Finally, in the case of Th > Tg then the spinodal transition may (but need not) interfere 
with our determination of Th- This is a significant ambiguity that we cannot resolve in the 
present calculations but the reader should be aware of its existence. 

We finish by listing some of the reasons that motivate our study. First it gives nonpertur- 
bative information about the free energy as a function of the Polyakov loop. While the value 
of Tc tells us when the ordered and disordered minima have the same free energy, the func- 
tion (Tofr(T) indicates how the curvature of the free energy at the confining vacuum changes 
with T and its vanishing may indicate the point at which the confining vacuum becomes 
completely unstable. This information can serve to constrain the form of the potential 



in effective models for the Polyakov loops (like those in for example). Second, this 
study is a first attempt to investigate the validity of nonzero temperature mean-field theory 
techniques of large- A^ lattice gauge theories in the continuum. A related study [li^ in the 
strong-coupling limit saw that mean-field theory predictions at finite temperature are simply 
incorrect, which is consistent with the fact that the large- A^ limit of these strongly-coupled 
fermionic systems is mapped to classical systems at finite temperature, whose universality 
class is not mean-field [ij]. It is of prime interest to know if the same happens in the contin- 
uum limit of the pure gauge theory, which we are much closer to. The Hagedorn temperature 
Th can also give an estimate of the central charge of a possible underlying string theory [20j , 
and finally it is interesting to know what is the limit of Th/T^ when A^ —>■ 00. This limit 
should be larger than 1, given that the deconfining phase transition has been shown to be 
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first order at 3 < < 8 jl5[- It is however possible that Th/T^ is close to 1 and perhaps 
decreasing with which opens up interesting new possibilities in the N = oo limit. 



II. LATTICE CALCULATION 

We work on a lattice with x Lt sites, where Lt is the lattice extent in the Euclidean 
time direction. The partition function is give by 

Z = J DUexp{-Sy^), (2.1) 

where S'w is the Wilson action given by 



1 - ^ReTr Up 

N 



(2.2) 



Here (3 = 2N'^/X, and A = g'^N is the 't Hooft coupling, kept fixed in the large- limit. 
P is a lattice plaquette index, and Up is the plaquette variable given by multiplying hnk 
variables along the circumference of a fundamental plaquette. Simulations are done using 
the Kennedy-Pendleton heat bath algorithm for the link updates, followed by five over- 
relaxations. We focus on the measurement of the correlations of Polyakov lines, which are 
taken every five sweeps. 

The Monte-carlo simulations were done using two different methods, which for conve- 
nience we refer to as method "A" and method "B" . In both, the calculation was initialized 
with a field configuration at the lowest value of (3 which was preceded by at least 3000 ther- 
malization sweeps from a cold start. Then, in method B, to reduce thermalization effects, 
each value of f3 was simulated beginning from the previous one. In contrast, in method A, 
more than one simulation began from the same (3. This serves to eliminate the dependence 
between different simulations^. We choose to use method A and method B for = 8 and 
A^ = 12 respectively. In the case of A^ = 10 we used both methods, and could compare 
between them. In Tables HlHTj where we present results obtained with method A, we give 
the initial configurations for each value of j3, along with the statistics of the study. The 
number of thermalization sweeps for A^ = 8, 10 for each value of (3 was always at least 3000. 
As mentioned, the SU{12) calculations were done with method B, and less thermalization 
sweeps (at least 400) were needed for each value of p. For each data set, we calculate the 
correlation functions of the thermal lines with improve d op erators jl7l |. and use a variational 
technique to extract the lightest masses 13, 2^ 2^ 2^, the results for which are given in 



lattice units, am, in Tables Hi ll V| where the errors are evaluated by a jack-knife procedure. 

To check that the results of method A are properly thermalized we also calculate the 
masses by excluding an additional 3000 sweeps (in addition to the standard thermalization 
sweeps). This results in a new set of masses, which we refer to as amhot, and which we 
present in Tables HHTTl We find that in general the thermalization effects are not significant, 
and in almost all cases am^ot agrees with am to within one sigma. 

The physical scale a^/a listed in the Tables was fixed using the interpolation for the string 
tension given in in the case of SU{8). For A^ = 10, 12 we extrapolate the parameters 



^ For a finite amount of statistics these two methods can lead to different results. In particular, method 'A' 
should be more noisy. 
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of the scaling function of ^3]; Co,i,2, and Pq/N^, in l/N"^ from their values at = 6,8. In 
addition we measured the string tension for N = 10 at f3 = 68.80 on an 8^ lattice, and for 
N = 12 at P = 99.2, 100.0, on an 8^*, and a 10"^ lattice respectively. The results, together 
with the string tensions calculated with the scaling function, are given in Table |3 We 
find that assuming an error only in the measured string tension, then the measured and 
calculated string tensions deviate at most by 1.6a. Evaluating the extrapolation error will 
make the agreement better. 

Let us note here that the points used for the fits exclude values of (3 in which tunneling to 
the deconfined phase was observed. Tunnellings were identified by observing double peaks 
in the histogram of the expectation values of the Polyakov lines, and by identifying clear 
tunneling configurations (again using the Polyakov lines as a criterion). For other values of 
/3, we observe an increase in the fluctuations of the lines with T, but did not see any tunneling 
configurations. Correspondingly, the histograms widen with T, and for some values of /3 start 
to develop asymmetric tails. Nevertheless in all these cases the overall average of the bare 
Polyakov line was always lower than 3 x 10~^. The (3 values of the excluded deconfining 
simulations are p = 44.01,44.015 for SU{8), and p = 69.10,69.14,69.16, and 69.17 for 
SU{10) (for method A), as noted in Tables HHITl For the case of SU{12) the last point of 
P = 99.95 was still confining after 10, 000 sweeps. At the next value of P = 100.00 the 
12^ X 5 system already showed signs of instability. 

Finally, to check the effect of finite volume corrections, which potentially increase as the 
mass decreases towards Th, we perform several additional calculations of the correlation 
lengths. In the case of A^ = 12 we performed a single calculation on a 16^ x 5 lattice at 
the smallest volume with P = 99.95. The mass obtained after 5,000 thermalization sweeps 
and 10,000 measurement sweeps, is am = 0.187(14), and agrees very well with the result 
obtained on the 12^ x 5 lattice (both are presented in Table HV|) . For A^ = 10 in method B, 
all values of P were simulated on both a 12^ x 5 lattice, and a 14^ x 5 lattice. The results 
are presented in Table ITTTl and we find at most a 1.3a difference between them. Comparing 
with the situation close to the second order phase transition of the SU{2) group [l7[, we find 
that for N = 2, finite volume effects are much more important than for A^ = 12. This puts 
the results of this work, which were largely done for only one lattice volume, on a more solid 
footing. This is also consistent with standard theoretical arguments that predict smaller 
volume corrections for gauge theories with larger values of A^. 

If we define the effective string tension, aes{T), to be the coefficient of the leading linear 
part of the free energy of two distant fundamental sources, then aeg{T) = mT where m is 
the lightest mass that couples to timelike Polyakov loops. Therefore the ratio ac^/a will be 
given by 

Oeff ^ am{p) 
a U-{a^{P)f 

At each simulated ratio of P we estimate T/T^ using 



% ay/a{P)'' 



(2.4) 



where {ay/a)c is the value of a^/a at P = Pc for Lt = 5, i.e. 5a{Pc) = T~^. It is extracted 
from Tc/y/a, which is 0.5819(41) for N = 8. (This assumes that T^j \fa varies at most 
very weakly with a(/5), which is in fact what one observes J7].) The corresponding value for 
A^ = 10, 12 is found by extrapolating in 1/A^^ according to measured values of Td \faiLt = 5) 
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TABLE I: Statistics and results of the Monte-Carlo simulations for = 8, obtained with method 
A. The "D" denotes tunneling to deconfining configurations, which is not considered for the fits. 
am hot is the result of the data set that excludes the first 3000 measurement sweeps. 



N = 8 


P 


am 
(all sweeps) 


am hot 


a^/a 


Initial 


(no. of sweeps)/10'^ 


43.850 


0.361(17) 


0.362(17) 


0.3615 


Frozen 


20.0 


43.875 


0.336(22) 


0.362(22) 


0.3580 


43.850 


15.0 


43.900 


0.334(17) 


0.337(17) 


0.3547 


43.875 


13.0 


43.930 


0.272(19) 


0.270(19) 


0.3507 


43.875 


17.0 


43.950 


0.314(15) 


0.326(15) 


0.3481 


43.850 


24.0 


43.975 


0.286(17) 


0.283(17) 


0.3448 


43.950 


10.0 


43.980 


0.258(18) 


0.245(18) 


0.3442 


43.975 


13.0 


43.985 


0.304(14) 


0.294(14) 


0.3436 


43.975 


13.0 


43.995 


0.206(13) 


0.207(13) 


0.3423 


43.975 


14.0 


44.000 


0.264(17) 


0.260(17) 


0.3417 


43.850 


29.0 


44.01;D 






0.3404 


44.000 


13.0 


44.015;D 






0.3398 


44.000 


13.0 


44.020 


0.240(17) 


0.253(17) 


0.3392 


44.000 


15.0 


44.025 


0.228(17) 


0.201(17) 


0.3385 


44.000 


16.0 


44.033 


0.213(14) 


0.220(14) 


0.3376 


44.000 


20.0 



for = 4,6,8 0. This gives 0.5758 for SU{10) and 0.5735 for SU{12), with an error of 
about 1%. 

In Figs. we give the effective string tensions as a function of temperature for the 
studied gauge groups and various data sets, as explained above. In view of Eq. ()1.4j) we 
choose to present (Tefr/a. As the relative error of am is roughly 10 times the one on Tcl\fo., 
we neglect the latter in our error estimate. The fit to the data was done according Eq. p.4p . 
either by fixing v = 0.6715, 0.5 for the 3D XY and mean-field universality classes, or in some 
cases by making u a free parameter. The fitting results are presented in Table IVII 

In Figs. lUHTTl we present confidence levels in the fit parameters, for the cases where the 
fit is good. In the case of two parameter fits, we present contours in the {A, Th/Tc) plane of 
the per degree of freedom levels, which correspond to confidence levels of 68.27%, 90%, 
and 99%. These confidence levels are then reflected in the error estimates we give in the 
text. 

When we treat i/ as a free fit parameter as well, we present two dimensional projections 
(e.g. in v— Th/Tc space) of a volume in the parameter space of (Th/Tc, A, u) that corresponds 
to the a confidence level of 68.27% and lower. In this case we do not quote in the text any 
error estimates together with the central values. 
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TABLE II: Results and statistics of the Monte-Carlo simulations for = 10 obtained with method 
A. The "D" denotes tunneling to deconfining configurations, amhot are the results of the data set 
that excludes the first 3000 measurement sweeps. 



= 10 





am 


amhot 




Initial 


(no. of 




(all sweeps) 






config. 


sweeps)/10^ 


68.5000 


0.571(41) 


0.609(34) 


0.3941 


frozen 


17.0 


68.5520 


0.524(28) 


0.514(29) 


0.3888 


68.50 


22.0 


68.6000 


0.545(11) 


0.544(11) 


0.3839 


68.50 


23.5 


68.6553 


0.502(13) 


0.500(14) 


0.3785 


68.60 


20.5 


68.7000 


0.481(36) 


0.441(35) 


0.3742 


68.50 


23.5 


68.7500 


0.452(12) 


0.453(14) 


0.3695 


68.50 


23.0 


68.8000 


0.391(17) 


0.388(19) 


0.3649 


68.50 


26.0 


68.8500 


0.389(14) 


0.378(22) 


0.3603 


68.50 


24.0 


68.9000 


0.356(15) 


0.345(16) 


0.3559 


68.50 


23.0 


68.9500 


0.287(15) 


0.289(16) 


0.3516 


68.50 


23.5 


69.0000 


0.295(17) 


0.286(17) 


0.3473 


68.50 


25.0 


69.0100 


0.295(13) 


0.294(14) 


0.3465 


69.00 


23.0 


69.0200 


0.315(15) 


0.312(15) 


0.3457 


69.00 


22.0 


69.0412 


0.277(15) 


0.292(16) 


0.3439 


69.00 


22.0 


69.0500 


0.271(14) 


0.273(16) 


0.3432 


69.00 


21.0 


69.0700 


0.251(18) 


0.242(19) 


0.3416 


69.00 


22.0 


69.0816 


0.262(16) 


0.258(17) 


0.3406 


69.00 


22.0 


69.1000 


_ 


_ 


0.3391 


68.50 


11.0 


69.1100 


0.253(11) 


0.248(11) 


0.3383 


69.10 


22.0 


69.1200 


0.233(15) 


0.226(16) 


0.3375 


69.10 


22.0 


69.1300 


0.223(15) 


0.231(14) 


0.3367 


69.10 


22.0 


69.1400 






0.3359 


69.10 


12.0 


69.1500 


0.218(12) 


0.228(13) 


0.3351 


69.10 


22.0 


69.1600 






0.3344 


69.15 


8.0 


69.1700 






0.3336 


69.15 


8.0 


69.1800 


0.222(14) 


0.225(16) 


0.3328 


69.15 


18.0 


69.1900 


0.218(12) 


0.228(13) 


0.3320 


69.15 


18.0 



III. SUMMARY AND DISCUSSION 

Gathering all the statistically reliable results of S'[/(10), we find that Th/T^ = 1.104(6) — 
1.114(15) when fitting with a 3D XY exponent. When u was treated as a free parameter 
fit we find central values of = 1 and Th/Tc = 1.160 — 1.172. Finally when a mean-field 
exponent is used, one find from data obtained with method B that Th/Tc = 1.087(11). We 
also find that the extrapolation with the mean- field exponent has a higher x^, both in the 
case of method A, and in the case of method B on a 14^5 lattice. In particular, when we 
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TABLE III: Masses in lattice units and statistics of the Monte-Carlo simulations for = 10 
obtained with method B on 14^ x 5, and 12^ x 5 lattices. 



iV = 10 




12^ X 5 


(no. of sweeps)/ 10^ 


14^ X 5 


(no. of sweeps) /lO'^ 


ay/a 


68.50 


0.599(30) 


20 


0.584(35) 


20 


0.3941 


68.70 


0.396(17) 


22 


0.428(29) 


22 


0.3742 


68.85 


0.371(27) 


20 


0.371(18) 


20 


0.3603 


68.95 


0.316(16) 




0.308(11) 


a 


0.3516 


69.07 


0.246(12) 


U 


0.270(15) 


(( 


0.3416 


69.13 


0.244(17) 


U 


0.241(14) 


18 


0.3367 


69.17 


0.220(10) 




0.203(15) 


16 


0.3336 


69.20 


0.154(11) 


22 


0.200(21) 


8 


0.3313 


69.23 


0.183(19) 


15 






0.3290 



TABLE IV: Results and statistics of the Monte-Carlo simulations for A^ = 12 obtained with method 
B. All masses are calculated on a 12'^ x 5 lattice, except for the last row which is for 16'^ x 5 lattice. 
The number of thermalization sweeps was at least 400 between two successive values of /3, while it 
was 800 for the first calculation of /3 = 99.00. 



N = 12 


/? 


am 


a^/a 


(no. of sweeps)/10'^ 


99.00 


0.573(25) 


0.3866 


5 


99.20 


0.472(24) 


0.3729 


5 


99.40 


0.396(18) 


0.3601 


5 


99.60 


0.333(17) 


0.3479 


10 


99.80 


0.259(13) 


0.3365 


10 


99.90 


0.203(13) 


0.3310 


10 


99.95 


0.190(13) 


0.3284 


10 


99.95 


0.187(14) 


0.3284 


10 


analyze both data 


sets together, we 


find that the best fit has a x^/dof 


= 1.93,1.36 (with 



do/ = 29) when fitting with a mean-field, and a 3D XY exponent respectively. It is also 
interesting that the most sizable contribution to the in this case comes from the data 
point at /? = 68.95 (also line no. 10 in Table ITT| . Since no reasonable fit will go through this 



TABLE V: String tensions for 5C/(10) and SU {12). 









Measurement 


Results of scaling function (see text) 


N 


= w,p = 


68.80 


0.3667(80) 


0.3649 


N 


= 12,(3 = 


99.20 


0.3770(26) 


0.3729 


N -- 


= 12,/3 = 


100.00 


0.3243(23) 


0.3257 
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TABLE VI: Results of the fits to Eq. ()1.4() . The values of Th/Tc^ and A are presented for all fits. 
In the cases where a fit with as a free parameter was made, then the resulting v is presented in 
the second column. We present the coverage probability contours in the A — Th/T^ plane, when 
appropriate, in Figs. IH tilll 



A T 

N 


Universality class 


Tc 


Th 


A 


central charge 


dof 


1 J? 

dof 


8 


3D XY 


1.093 


0.636 


2.190 


1.178 


2.84 


11 


12^ X 5 


MF 


1.066 


0.620 


1.728 


1.237 


2.79 


11 


(method A) 


Free, v = I 


1.145 


0.666 


3.055 


1.073 


3.19 


10 


10 


3D XY 


1.104 


0.636 


2.328 


1.236 


1.44 


21 


12^ X 5 


MF 


1.078 


0.621 


1.862 


1.178 


2.07 


21 


(method A) 


Free, = 1 


1.160 


0.668 


3.132 


1.067 


1.16 


20 


10 


3D XY 


1.108 


0.624 


2.112 


1.169 


3.29 


7 


12^ X 5 


MF 


1.083 


0.638 


1.681 


1.223 


3.56 


7 


(method B) 


Free, u = 0.79 


1.127 


0.649 


2.394 


1.131 


3.94 


6 


10 


3D XY 


1.114 


0.642 


2.102 


1.156 


0.68 


6 


14^ X 5 


MF 


1.087 


0.626 


1.683 


1.215 


1.00 


6 


(method B) 


Free, = 1 


1.172 


0.675 


2.815 


1.045 


0.60 


5 


12 


3D XY 


1.116 


0.640 


2.337 


1.162 


0.30 


5 


12^ X 5 


MF 


1.092 


0.626 


1.858 


1.215 


0.50 


5 


(method B) 


Free, u = 0.69 


1.119 


0.642 


2.387 


1.156 


0.37 


4 



point (see Fig. E}, it is quite conceivable that at this value of (3 we have a strong statistical 
fluctuation. Ignoring this point gives essentially the same fitting parameters, but makes 
X^/dof = 1.54,1.04 for the mean-field and 3D XY exponents respectively. This suggests 
that the 3D XY exponent is preferable (although we still cannot exclude the mean-field 
exponent possibility completely). This preference is also be seen in Figs. |71 and IHl where 
we give the projections in the z/ — Th/Tc plane, of a volume in the {u, A, Th/Tc) space, that 
corresponds to a confidence level of < 68.27%. There we see that the point u = 0.5 is either 
outside, or at the edge of this volume. 

For = 12 the fits are better, and we find that Th/Tc = 1.092(6) — 1.116(9) for mean- 
field and 3D XY exponents. Both fits have a low x^, and we again cannot rule out either. 
Here again the mean-field fit has a higher than that of the 3D XY fit. In this case, 
however, in view of the good values, the preference towards a 3D XY exponent is weaker 
than for SU (10). Nonetheless it is interesting that when we perform a fit with the exponent 
z/ as a free parameter we find Th/Tc = 1.119, and z/ = 0.69 which is closer to the 3D XY 
exponent z/ = 0.6715(3) than to the mean-field value of u = 0.5. Looking at Fig. ITD we find 
that the preference of our data for u = 0.6715 is indeed weaker here, as both exponents have 
a similar position in the parameter space, with respect to the shaded area. 

The limited statistics prevents us from making statements about the behavior of Th as 
a function of A^. This is unfortunate, since it is of interest to know how far is Th/Tc from 
1 at = oo. Nevertheless we obtain fitted values of Th/ \fo ~ 0.62 — 0.68, which is lower 
than Tc/ \fo ~ 0.7 of SU{2) ji^] where the phase transition is second order, and therefore 
may be Hagedorn, Tc = Th, or provides a lower bound on Th- To emphasize this point we 
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FIG. 1: Effective string tension for SU{8) obtained in method A in units of the zero temperature 
string tension. 



plot the o'cs/o' for = 10, 12 that we obtained with method B here. For guidance we also 
calculate the masses from Polyakov loops for SU{2) on a L^5 lattice with L = 16, 20, and 
24. The results are presented in Table IVIH and also show the expected increase in finite 
volume corrections as the temperature approaches Tc. For each value of T j \fo we choose to 
present in Fig. the corresponding value of o^'g^ja calculated on the largest volume there, 
together with an extra point for SU{2) at T = with (Jcfr = 0. For the physical scale we 
again use the interpolation fit in We find that for = 10, 12 the results seem to be 
close, but not to follow the SU (2) results. 

Finally we find interesting the fact that the Nambu-Goto action, gives rises to a Hagedorn 
behavior at similar temperatures given by 

T,^^/v^ ~ 0.691/v^, (3.1) 

where c is the central charge ji^l and equals unity in the usual Nambu-Goto model. Applying 
this formula to our values of Th we give the central charge values listed in Table IVII 

A determination of the proper universality class is an important issue, and may teach 
us how the scaling region behaves with N (if at Th the system behaves like in a second 
order transition). As mentioned in the introduction, this question was studied numerically 
in [l^ for the strongly-coupled gauge theory (with quarks included), and analytically in 
|l^ where the scaling region for chiral restoration was seen not to change with N. In our 
context the similar question can be approached for deconfinement in the continuum of the 
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FIG. 2: Effective string tension for SU{10) obtained in method A in units of the zero temperature 
string tension. 



pure gauge theory. Unfortunately, despite the mild preference towards the 3D XY model, 
discussed above, we currently cannot rule out unambiguously any of the universality classes. 
To distinguish which one is actually correct one must approach Th closer, and increase the 
statistics. In fact, as discussed in the introduction, if the Hagedorn transition is second 
order, we expect that the critical region shrinks with increasing A^, and that only when 
{Th — T) < one will see the nontrivial critical behavior of a 3DXY model. This 

is however hard because tunneling configurations become more probable, and finite volume 
effects (although relatively small at larger values of A^) become more important. Nonetheless 
when we perform a fit with the critical exponent u G [0, 1] as a free parameter we find that 
for SU (12), the best fit result for u is closer to the 3D XY exponent than to the mean-field 
value. 

We believe that a more thorough investigation (with larger statistics) would render the 
understanding of the proper universality class, and indeed all the other issues discussed 
above, much clearer, and the large- limit of the gauge theory better understood. However, 
considering its current numerical cost we postpone it to future studies. A different route to 
approach this question is to study variants of the pure gauge theory, such as adding scalar 
fields. Depending on the couplings added, the phase transition can become second order 
[iij l , and one can study the adequacy of large- iV mean-field techniques close to second order 
phase transition in a thermodynamically stable phase. 
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FIG. 3: Effective string tension for 5^7(10) obtained in method B on a 12'^ x 5 lattice, in units of 
the zero temperature string tension. 
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FIG. 6: Confidence levels for the fits of SU (10) on a 12^ x 5 obtained with method A for the 3DXY 
exponent. 
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FIG. 7: Two dimensional projections of a isosurface in the three-dimensional space of (v, Tfj/Tc, A) 
which represents a confidence level of 68.27% for the fits of SU{10) obtained with method A on a 
12^ X 5 lattice, with a free exponent. 
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FIG. 8: Confidence levels for tlic fits of SU{10) on a 14'^ x 5 obtained with method B. The plot on 
the right(left) is for the 3DXY(MF) exponent. 




FIG. 9: Two dimensional projections of a isosurface in the three-dimensional space of (z^, Th/Tc, A) 
which represents a confidence level of 68.27% for the fits of SU (10) obtained with method B on a 
14^ X 5 lattice, with a free exponent. 
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FIG. 10: Confidence levels for the fits of SU{12) obtained with method B, with mean-field theory 
and 3D XY exponents. The plot on the right is for the 3D XY exponent, and on the left with a 
mean-field theory exponent. 
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FIG. 11: Two dimensional projections of a isosurface in the three-dimensional space of (i^, Th /Tc, A) 
which represents a confidence level of 68.27% for the fits of SU{12) obtained with method B, with 
a free exponent. 



19 



0.8 



0.6 



" 0.4 



















T 1 

f * 










• 
■ 


JN=z 
N=12 


















^ N= 10, method B, 14 x5 




















T 1 


1 














^ 






















i: 




















































— • — 



0.2 



"0:5 



0.55 



0.6 



0.65 



T/a 



1/2 



FIG. 12: Effective string tension for N = 2, 10, 12 in units of the zero temperature string tension. 



TABLE VII: Masses from Polyakov loops on L^5 lattices for N = 2. 

loop mass: SU(2) 



(3 


16^5 


20^5 


24^5 


a^/a 


2.28 


0.472(11) 


0.476(13) 


0.460(26) 


0.3871 


2.29 


0.439(9) 


0.435(8) 




0.3754 


2.295 


0.417(22) 




0.383(22) 


0.3639 


2.30 


0.373(11) 


0.390(12) 




0.3527 


2.31 


0.359(8) 


0.335(8) 


0.368(16) 


0.3417 


2.32 


0.303(8) 


0.299(9) 


0.299(7) 


0.3400 


2.3215 


0.290(19) 




0.288(22) 


0.3309 


2.33 


0.274(8) 


0.252(8) 


0.245(8) 


0.3256 


2.335 


0.252(11) 




0.222(10) 


0.3204 


2.34 


0.217(8) 


0.205(7) 


0.196(7) 


0.3101 


2.35 


0.161(6) 


0.158(7) 


0.158(7) 


0.2891 
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